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DIRECTIONAL THERMAL-RADIATIVE PROPERTIES OF CONICAL CAVITIES 
by Leslie G. Polgar and John R. Howell 
Lewis Research Center 

SUMMARY 

A Monte Carlo technique is applied to the problem of finding the apparent directional 
thermal-radiative reflectance, absorptance, and emittance of a right circular conical 
cavity with a diffusely reflecting internal surface. Parameters examined are the cone 
angle, the cone surface absorptivity, and the angle of incident radiation. Small cone 
angles are found to produce a strong reflection back to the direction of incident radiation, 
while the results for large cone angles approach those for a flat plate. The Monte Carlo 
method is found to give a useful tool for solving problems of this type. The directional 
emittance results may be useful in the field of optical pyrometry, where the effect of 
viewing angle on the apparent temperature of a cavity is of interest. Some correlation to 
the reflection properties of the Moon are noted. 




INTRODUCTION 

The directional radiative properties of surfaces are nearly always neglected when 
calculations involving radiative -energy transfer are made. This is due both to the lack 
of experimental data on directional effects for real surfaces and to the difficulty of in- 
cluding such data in routine calculations. The reflection laws usually assumed are for 
diffuse (cosine law) or specular (mirrorlike) reflections, or some combination of these, 
as suggested in reference 1. The former mode assumes the reflections to be independent 
of the angle of incident radiation, while the latter assumes the simplest possible relation 
between angle of incidence and reflection, that is, the angle of incidence equal to minus 
the angle of reflection. Both cause a great reduction in the complexity of calculation. 

Most common surfaces do not reflect in either a specular or a diffuse manner. The 
Moon is one common example (ref. 2) where the reflectance* of visible radiation is 

*In this report, the ending "ance" (i. e. , reflectance) refers to apparent surface 
properties, whereas the "ivity" ending (reflectivity) refers to a physical property of the 
material. 



strongest back in the direction. of the incident radiation. More discussion of this is 
given in the DISCUSSION OF RESULTS. It has also been shown (refs. 3 and 4) that cavi- 
ties of various shapes can greatly alter the directional reflective and emissive properties 
of surfaces and can reflect energy back in the direction of incident radiation. 

Two-dimensional shapes can be analyzed by standard methods involving formulation 
of the energy exchange equations in terms of "shape” or "exchange" factors available in 
the literature for many simple geometries. Some shape factors for three-dimensional 
bodies can be found in the literature, but the complexity of formulation of radiative ex- 
change problems becomes great. For example, Sparrow and Jonsson (ref. 5) treated the 
relatively simple case of finding the total hemispherical emittance of a conical cavity as- 
sumed to be either at constant surface temperature or having uniform surface energy 
flux. 

In this report, a Monte Carlo method is used as an alternate approach to the usual 
shape factor formulation in finding the directional properties of a conical cavity, with a 
view to developing some insight into surface anomalies and developing a useful technique 
for solving problems of this type. The results are useful in determining the accuracy of 
pyrometrical measurements of the surface temperature when this temperature is meas- 
ured by viewing a cavity in the surface. 


SYMBOLS 


A 

B 

C 

G 

1 

P 

R 


area of element on hemisphere over reflecting surface 
cos £’/cos 9 

cos sin rj’ cos \{/'/ cos 9 
cot 9 / cot £ 

path length of bundle between reflections 

probability of incident energy bundle striking given illuminated area element 
on internal surface of cone 

number chosen at random from evenly distributed set of numbers in range 
0 to 1 


r distance along element of cone 

X, Y, Z Cartesian coordinate system defined in fig. 7(a), Z being on or parallel to 
the cone axis and X being the line of intersection of the cone base and the 
plane perpendicular to the base and parallel to the incident radiation 

X'^’jZ’ coordinate system defined in figs. 2 and 7(a) 
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a a apparent absorptance of cavity 

a g surface absorptivity 

y angle between lines in base plane parallel to projections of bundle path and 

incident energy, fig. 1 

6 incremental value 

C angle between cone axis and line through cone axis parallel to incident radia- 

tion, angle of incidence 

rj angle between normal to cone wall at point of bundle path intersection and 

bundle path 

6 cone half-angle 

M local cone radius at point of bundle reflection 

p directional reflectance of conical cavity 

<p circumferential angle around cone axis, measured counterclockwise from 

X-axis to p 

\p angle from X’ -axis to projection of bundle path in plane tangent to cone wall 

at point of bundle intersection 

£2 variable 

Subscripts: 

i, j angular increment indices 

max maximum 

min minimum 

p projected normal to direction of incident radiation 

£,£’, 7 ’ function of directions and y’, respectively 

Superscript: 

1 values after reflection 

ANALYSIS 

Incident parallel radiation, assumed to consist of discrete bundles of energy, is im- 
pinging at a given angle on a right circular conical cavity of known surface absorptivity. 
The base radius of the cone is used to nondimens ionalize all lengths in the analysis. The 
walls of the cone are assumed to be gray and diffusely reflecting; polarization effects 
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are not considered. It is assumed 
that no energy is emitted from the 
cone walls (i. e. , the cone is cold 
compared to the source of incident 
radiation). 

A Monte Carlo method is used to 
find the apparent absorptance and the 
distribution of reflected energy from 
the conical cavity. This method con- 
sists of following each discrete bundle 
of incident energy through its probable 
path in the region of the cone, taking 
into account the diffuse reflections 
within the cone and the absorption of 
energy bundles at the gray internal 
surface. Those bundles which leave 
the cone are tallied in the angular 
increment Ayl, A£t on an imaginary 
hemisphere over the mouth of the 
cone (fig. 1) from which they leave. 

The initial point (t±, <p) where the 
bundle strikes and the question of 
whether or not the bundle is absorbed 
or reflected at that point are both 
determined by a Monte Carlo tech- 
nique. If it is absorbed, this is 
tallied. If it is reflected, the direc- 


tion of reflection is specified by two angles, \p' and rj’ (fig. 2), that are chosen at ran- 
dom from the properly weighted set of angles of reflection. It is determined whether the 
reflected bundle exits from the cone or traverses the cavity and intersects the cone wall 
at another point. Exiting bundles are classified by their angular coordinates (y’ and £’). 
Bundles that strike the cone wall a second time are again checked to see if absorption 
occurs at this point. The internal reflection process is continued until the bundle is 
either absorbed or leaves the cavity. 

The apparent absorptance a a defined here as the fraction of energy incident 
from a given direction £ that is absorbed within the conical cavity, is then simply the 
fraction of incident bundles that is absorbed within the cavity. Defining the directional 
reflectance p r , as the fraction of energy incident from direction £ that is re- 
fleeted per unit solid angle around , y' allows determination of this value by tallying 
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the fraction of total incident bundles that leave in each unit solid angle. 


Point of Initial Bundle Intersection 

The illuminated area of an element on the conical surface projected normal to the 
direction of incident radiation is given by 

dA = cos r\ dA (1) 

where the area of an element of the conical surface is given by 

dA - [i esc 6 dfi. dtp (2) 

and cos t) is found in terms of known quantities by projecting a unit vector of incident 
radiation onto the normal to the cone wall at the point of intersection. This is done di- 
rectly and by components (appendix A, eq. (A7)) to give 

cos tj = sin 0 cos £ + cos Q sin £ cos tp (3) 

The probability of an incident energy bundle striking a given illuminated area ele- 
ment on the internal surface of the cone is the projected area of that element over the 
total projected illuminated area: 


P(dM dtp) = c ? s 71 ^ (4) 

/ cos rj dA 
J k 

This probability can be separated into the probabilities of a bundle striking an element 
dp and an element dtp because substitution of equations (2) and (3) into equation (4) in- 
dicates that ii and <p appear only as products and are therefore independent in the 
probability sense. The separation gives 


P(dM) = /^ maX P(dM d ^ ) 
v min 


( 5 ) 


and 
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Cone 

axis 

Z 


( 6 ) 



Incident 

radiation 


Figure 3. - Graphic definitions of M m | n , <p m in 
and <p maj( . The boundary of the area visible 
to the incident radiation is defined by /j. min 
which is the least /j. defined for a given cp, 
where < Pmin^‘P$° or °<'P<<Pmax- 


and a must be calculated. 

lUdX 


P(d <p) = P(d m dtp) 


mm 


When the angle of incident energy is less than the 
cone half-angle, that is, when £ < 6 , the entire inner 
surface is illuminated. Consequently, the local radius 
/i m^ at the minimum illuminated point for any value 
of cp is 


u . = 0 
'mm 


(7) 


and the minimum and maximum circumferential angles 

w and w . at which illumination occurs are 
r max r min 


^max “^min 71 


( 8 ) 


For £ > 0, part of the conical cavity is shadowed by 
the lip of the cone (fig. 3). For these cases, P min 
From the relation 



2 2 
juT . sin^w ) 
'min ■ j 


+ Vin COS <P = 


f 1 - >W cot 9 

tan(90 - C) 


(9) 


derived from relations for the depth at which the incident radiation strikes the wall at a 
given cp, and for a given angle of incidence (appendix A, eq. (All)) M min is given by 


M 


min 


G(G + cos cp) + (G cos 3L±J1 
G 2 + 2G cos cp + 1 


( 10 ) 


where G = cot 6 / cot The top sign is the correct choice as examination of any special 
case with P min ^ 1 will indicate. This choice leads to 


M 


min 


G 2 - 1 


G^ + 2G cos cp + 1 


(ID 


To find <? max (always equal to <P min by symmetry), P min is set equal to 1 in equa- 
tion (11). Then 
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Figure 4. - Determina- 
tion of distribution of 

Q. 


cos '"max = ^ < 12 > 

There remains only the problem of choosing values of /i. 
and <p such that the distributions given by equations (5) and 
(6), respectively, will be satisfied after enough values are 
chosen. 

For functions that can be easily integrated, a straight- 
forward procedure exists for doing this. The procedure con- 
sists of carrying out the integration of the probability P(ft) 



and choosing values of R in a random fashion to give values of 
the variable ft. The physical meaning of this procedure is demonstrated by figure 4. 
Because the integral of P(ft) over all ft must be unity, R can only take on values in the 
range 0 to 1. Equation (13) gives the probability of any element P(ft)dft lying in the 
range ^ min to ft. Consider choosing five equally spaced values of R. If in each equal 
increment AR an equal number of values of R is chosen, for example N, the value 
N/Aft in each increment will be proportional to the slope in that increment. But the 
slope in each increment must be directly related to the probability P(ft) corresponding 
to ft for that increment since dR/dft at ft is exactly P(ft). Choosing enough values 
of R in small enough intervals and making sure they come from a set of R evenly dis- 
tributed in the range 0 to 1 guarantee accurate derivatives. Therefore, the distribution 
of ft calculated will adequately describe the original distribution P(ft). 

Equation (13) can be applied for any integrable probability distribution. For j u and 
c p this process gives 



(14) 


and 
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r<p 

Jrrt 


R„ = 

<P 


‘'cp 


(cos £ + cot 6 sin £ cos cp)dcp 


min 


f max (cos £ + cot 6 sin £ cos cp)dcp 
<P r 


min 


(15) 


respectively. Solving for ju and <p gives 


^ = ” ^min) + ^minl 


1/2 


(16) 


and 


<P = - 1 )(<P max + c ot 6 tan ? sin <P max ) - cot 6 tan £ sin cp 


(17) 


If two numbers are chosen at random from a set of numbers evenly distributed between 
0 and 1, and they are designated as and R , an initial point cp at which an incident 
energy bundle strikes the interior surface of the cone can be chosen from equation (17). 
Then M min is found from equations (10) or (7) as the case may be, and /! is evaluated 
from equation (16). 

Equation (17) is transcendental in form, involving both cp and sin cp, and its solu- 
tion requires iterative means. Because of the rapidity with which this program as a 
whole is expected to run, no great pains are taken to minimize running time for this 
step, and the Newton-Raphson technique, modified to include third-order terms, is em- 
ployed (see appendix B). An iteration check is used, and if convergence is slow, an al- 
ternate regula-falsi method is used (ref. 6). 


Determination of Bundle Absorption 

When the original point (i u, cp) at which the bundle strikes the wall is determined, the 
probability of bundle absorption at that point is checked by comparing another randomly 
chosen number R to the surface absorptivity a of the cone. If 

CL S 

R a < 0! s (18) 

the bundle is tallied as absorbed, and a new incident bundle history begun. Otherwise, 
reflection must occur. 
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Angles of Reflection From Cone Surface 



Figure 5. - Reflection from cone 
mouth. 


The angles of reflection from a diffuse wall are determined 
by the relations 


\f/ = 27tR^ (19) 

cos r\ = VS” (20) 

as derived in reference 7 for reflection from a diffuse surface. 


Determination of Bundle Destination 

If the bundle is reflected from the position \i, cp on the cone 
surface at angles tj* and ip ' , the next step is to determine its 
destination. The angle to the cone axis at which the bundle is 
reflected gives one convenient check (fig. 5), for the bundle must 
leave the cone mouth if 


cos < cos 6 


(21a) 


where the equation for cos C' is found by resolving £ into its components along X', 

Y' , and Z' , projecting the components along Z, and then equating the sum of these com- 
ponents with the direct projection of £ along Z. This procedure gives 


cos = sin 6 cos rj' + cos 6 sin tj’ cos \p' 


(21b) 


If > 6 , the bundle may or may not leave the mouth, however, and a further check is 
needed. 

If the bundle may hit the cone wall, this fact can be determined by finding the point 
(M',^') on the wall or its extension where it will strike. The determination is discussed 
in the next section. If 


M' < 1 


( 22 ) 


the bundle will strike the cone surface. 
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Determination of Point of Intersection With Cone Surface 


The value of the point of intersection (p 1 , cp') must be determined in terms of the last 
point of reflection (ju, cp) and the angles of reflection from that point (r)' 

In appendix A, it is shown that 

cos(Zr,r') = cos 2 0 + sin 2 0 cos (cp' - (p) (23) 

When equation (23) is substituted into the law of cosines for certain imaginary triangles 
within the cone, it is shown further in appendix A that two independent equations are ob- 
tained: 


jf 2 = r 2 + r’ 2 - 2rr'[cos 2 e + sin 2 0 cos (cp' - <pj] 
i 2 sin 2 £' = sin 2 ^’ 2 + r 2 - 2rr' cos(cp' - cp^ 


(24) 

(25) 


A third relation can also be found: 


2 2 2 

r’ = r + SL + 2r£ sin rf cos \p' 


(26) 


When equation (21b) and equations (24) to (26) are used, cp ' is obtained as 


cp' = cp + cos 


r - r’ cos 2 9 ± ( r< - r ) sin cos V cos e 
-1 cos C’ 


r* sin 2 0 


(27) 


with choice of the top sign being dictated by examination of the special case of cp' - cp = it, 
6 = 30°, rj* = 29, and = (tt/2) + 6. 

Considerable algebraic manipulation gives 


r , _ r 1 - C ± (B* - C) 
(1 - B 2 ) 


(28) 


where 


b - cos V 

cos 6 


10 



and 


C _ cos g* sin rf cos ij ' j' 
cos 9 

Examination of the special case of reflection in a direction normal to ju, q> shows that 
the plus sign must be chosen in equation (28). 

The value of fi' is of course found directly from r’ by the relation 


p.’ - r’ sin 9 (29) 

and substitution of the value of r' into equation (27) gives cp' . If p.', cp' is indeed on 
the cone surface, it is again determined whether absorption will occur. If the bundle is 
not absorbed, new angles of reflection are determined, and this process is continued un- 
til the bundle is either absorbed at the cone surface or reflected from the mouth of the 
cone. 

If the bundle is reflected from the mouth of the cone, the angles of reflection with 
respect to the cone axis (y* and £’) must be determined (fig. 1, p. 4). 


Determination of Angles of Reflection From Cone 

The angle of reflection can be computed by the equality given in equation (21). 
The angle y’ can be computed from the value of cp' at the last point of reflection 
(H' ,cp') and the angles of reflection from that point (rj’ ,i|/ T ). This relation is derived in 
appendix A as 


i t „„„-l cos 5* sin 9 - cos n’ 

Y - <p' - cos ^ (30) 

sin £’ cos 9 

Care is taken in using this relation to avoid the problems inherent in allowing the digital 
computer to choose the principal values in evaluating trigonometric functions. 

The angles at which the bundle leaves are then examined to find which increment 
(Ay* , A^’) it passes through on the surface of the imaginary hemisphere over the cone, 
and is tallied in that increment. After a sufficient number of bundle histories are fol- 
lowed, the values of the directional properties are computed. 

The computer program is written in the Fortran IV language and run on the IBM 7094 
digital computer. The program is relatively straightforward, as can be deduced from 
the flow chart (fig. 6, p. 12). 
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Figure 6. - Computer flow chart. 
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(b) Effect of cone half-angle; cone surface absorptivity, 0. 25. 
Figure 7. - Apparent absorptance. 


Running time is primarily a function of 
the number of histories followed and the cone 
half-angle 9. The half-angle affected the 
number of internal reflections. For 6 of 
1/2° and any given angle of incidence and sur- 
face absorptivity, the running time for 
50 000 histories was about 9 minutes. For 
the same conditions but with 9 of 89. 9°, 
running time was 3. 5 minutes. Surface ab- 
sorptivity also has some effect; for an almost 
black surface, most bundles are absorbed on 
the first encounter with a wall, and few re- 
flections need be followed. 

DISCUSSION OF RESULTS 

The apparent directional absorptance of 
a conical cavity is shown in figure 7. As has 
been derived in reference 4, the directional 
absorptance is equal to the directional emit- 
tance for any isothermal cavity. This is 
true for monochromatic radiation or, for a 
cavity with a gray surface, for the total di- 
rectional absorptance. 

The apparent absorptance for a cone of 
any angle is always greater than the cavity 
surface absorptivity because of multiple re- 
flections within the cone. The maximum ap- 
parent absorptances o' r occur for the 

a > s 

smallest angles of incident radiation £ . For 
large angles of incidence, a? a ^ diminishes 
because the fraction of incident energy under- 
going multiple reflections is reduced. Only 
that part of the cavity wall near the lip is 
illuminated by the incoming radiation. 

Figure 8 depicts the directional reflec- 
tance of a cone with a 15° half -angle for 
radiation incident at 60° to the cone axis. 
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Angle of reflection in plane parallel 
to incident radiation, 

V, 

deg 


Surface 

absorptivity, 

“s 

0 



The results in figure 8 are shown only in the plane containing the cone axis and the di- 
rection of incident radiation, and the effect of surface absorptivity is demonstrated. 

In figure 9, the directional reflectance of perfectly reflecting cones with various 
cone half-angles between 1/2° and 89. 9° is presented for four angles of incident radia- 
tion. Results are given in two planes, each perpendicular to the base of the conical cav- 
ity and containing the cone axis. One plane is parallel to the direction of incident radia- 
tion, and the other is perpendicular to that plane. The most striking result is the pres- 
ence of strong reflections occurring in the direction of the incident radiation, especially 
at the smaller cone angles. A secondary peak occurs opposite the strong reflection be- 
cause of the internal reflections within the cone. The secondary peaks should decrease 
rapidly in magnitude with increasing surface absorptivity. As the cone angle is in- 
creased, the results approach the curve for a diffuse flat plate as expected. 

Figure 10 (p. 17) presents the circumferential variation of reflectivity for various 
combinations of parameters, again showing the generally preferred reflection to the di- 
rection of incident radiation. 

Statistical variations in the Monte Carlo results are pictured in figure 11 (p. 18) for 
a representative case to indicate the probable accuracy of the results. Because the an- 
gular increments become quite small near the normal to the cavity base, the number of 
Monte Carlo energy bundles passing through these increments is small. The small num- 
ber leads to larger expected deviations near the normal. 

The interesting parallel between the reflectivity results of conical cavities and those 
for the surface of the Moon is noted in figure 12 (p. 18). The observed lunar values 
given in reference 2 are compared to the normalized reflectivity results computed by the 
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Angle of reflection in plane parallel Angle 0 f 

to incident radiation, incidence. 



(a) Cone half-angle, 1/2°. <b> Cone half-angle, 1/2°. 



Figure 9. - Directional reflectance of conical cavity. Cone surface absorptivity, 0. 





Angle of reflection in plane parallel 
to incident radiation, 
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V. 
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Figure 12. - Normalized directional reflectance of lunar surface from reference 2 compared to conical cavity re- 
sults. Cone half-angle, 15°; surface absorptivity, 0 . 5 . 



method of this report for a conical cavity with a half-angle of 15° and a surface absorptiv- 
ity of 0. 5. The similarity of these results may indicate some facts about the fine struc- 
ture of the lunar surface. First, small surfaces near the normal might be expected to 
exist because of the relatively small cone angle which provides the similarity. Second, 
the lunar craters of small size which were observed by photographs from the Ranger 
lunar probes probably do not contribute greatly to these reflectivity effects. This is be- 
cause the crater half-angles are much greater than 15°, more often approaching 60°. 

The reflectivity effects must therefore be based on much smaller structure. 


CONCLUSIONS 

The apparent absorptance and directional reflectance of conical cavities with diffusely 
reflecting internal surfaces are determined by a Monte Carlo procedure. The Monte 
Carlo technique is found to be a useful tool for solving problems of this type and appears 
to offer promise, especially for complex geometries. 

Analytical results indicate that a strong preferential reflection in the direction of the 
incident beam of radiation is present and is most apparent at small cone angles. Because 
the only limitation on the size of cavities causing these effects is that they be large com- 
pared to the wavelength of incident radiation, it is possible that surface microstructure 
could strongly affect the reflectance of engineering surfaces. Measurements made on 
polished laboratory samples will probably not predict the directional behavior of less 
highly finished surfaces. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 14, 1965. 
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APPENDIX A 


DERIVATION OF BUNDLE PATH EQUATIONS 
Relations Between Incident Angle 


In figure 13, the interrelations between angles for an incident bundle are shown. To 
derive an expression for the angle r? in terms of the angles £, 0 , and <p, the length 
C - S' can be divided as 


(C - S') = {C - S) + (# - S') 


Then 


{0 - S') = {€ - J ) cos t\ 


(Al) 


(A2) 


and 


( 0 _ j) - (* - y) - (* - *0 - (* - -^) sin j cos <P 
cos 0 cos 8 cos 8 


(A3) 


Also, 


where 


(S - S') = S)sin 6 = [(j* - V) - (S - 9')] sin 8 


(A4) 


(j*- <g) = {0 - »cos £ (A 5) 

and 

( & -<?) = ( 0 - S )tan 0 = (0 - J )sin £ cos <p tan 6 (A6) 

Substituting (A2) to (A6) into (Al) gives 

cos rj _sin_£_cosjp + cog ^ g ^ n g _ S j n ^ cog ^ g ^ n g j. an g _ cog ^ S j n g + S i n £ cos q cos ^ 
cos 8 

(A7) 
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Derivation of p m j n Relations 


The minimum point on the interior surface of the cone at which a bundle may strike 
for a given value of the angle cp can be derived by referring to figure 14. The line 
(s/ - <g ) has a length given by 


(j/- <r) = (3> - <r)cot e = (i - M min )cot e 


or, alternately 


(j* - V) = {9 - tf)tan(90 - 0 


(A8) 


(A9) 


where 


(9 - <g) = (V -f) + (f- 9) = ^x min COS + (1 - M^ in sin 2 (p) 1/2 (A10) 


Substituting (A 10) into (A9), and equating the result to (A8) gives 


Cone 

axis 

Z 


(1 - M min )cot 0 = tan(90 - £) 




v 


Figure 15. - Derivation of bundle path relations. 


T 2 2 V2] 

x l/'min cos V + < X - ^min sln J ( A11 > 

Derivation of Bundle Path Relations 

The bundle path relations may be found from 
figure 15 by writing the law of cosines for triangles 
(s/09), (90 g), and (s/9 c g), to give, respectively 

£ 2 = r 2 + r 2 - 2rr' cos(Z r, r’) (A12) 

(9 - g ) 2 = I 2 sin 2 Z,' = r 2 sin 2 0 + r’ 2 sin 2 0 

- 2rr' sin 2 9 cos (cp* - cp) (A13) 

£ 2 = - <^) 2 + (r cos 6 - r f cos 0) 2 (A14) 

Combining to eliminate £ 2 and (9 - g) gives 
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Figure 16. - Auxiliary coordinates along cone wall 



Figure 17. - Choice of sign in 
equation (27). 


2 2 

cos(Zr,r') = cos 6 + sin' 6 cos (cp' - cp) 

Also, using the coordinates shown in figure 16, the law of cosines gives 

r’^ = r^ + J - 2xH cos(Z It, r) 


but 


(A15) 


(A16) 


cos( Z Jfl, r) = - ^ sin cos ^ 
£ 


and substituting (A17) into (A16) gives 

2 2 2 

r' - x + SL + 2r i sin t]' cos \p' 


(A 17) 


(A 18) 


Correct Choice of Sign for Equation (27) 

Figure 17 shows the special case of cone angle 6 = 30°, incident angle 
£ = ff/2 + 6 = 120°, and the reflected bundle passing through the cone axis. It follows 
that 7]’ = 26 — 60°, if/' = 77, and (cp 1 - cp) = ir. Also, the point of intersection on the cone 
wall, r' , is related to r by 


r’ = r sin 6 = — 
2 


(A19) 
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Cone 

axis 

Z 



Substituting into equation (27) gives 


r« CO s 2 e 1 fr' ~ r ) sin V' cos cos 6 


cos (<p’ - <p)= -1 = 


cos 


r' sin 2 0 



(_r')^/3\/3 (_!) 

2r' - 2r * ± 

2 2 

4 

fi) 


- \ 2/ . 


4 


(A20) 


or 


1 

4 



(A21) 


from which it follows that the plus sign must be chosen. 


Determination of Bundle Exit Direction 


From figure 18, it can be seen that the direct projection of the length (.sZ-tf) = 1 
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directly onto (jtf - <S) is simply cos rj. Projection of (st - <Jf) onto {s/ - <g) by components 
gives 

(sf - )cos t) = (s/ - <i?)[cos sin 8 + (-cos(<p f - y’)sin cos 0)] (A22) 

or 


y« 


fpt _ cos -l cos C’ sin 8 - cos r\ 
sin cos 8 


(A23) 
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APPENDIX B 


SOLUTION TECHNIQUE FOR TRANSCENDENTAL EQUATION 

The transcendental equation in cp (eq. (27)) was solved by using a modified Newton- 
Raphson technique. The standard Newton-Raphson expression is 


<P = <P o 


% 0 ) 


(Bl) 


where cp q is the value of cp from the previous iteration. When a third-order term is 
desired, the form is 


% 0 ) 

(p ' (p ° = 'pu ~) {1 + 5) 


(B2) 


The Taylor series expansion of cp when f (cp) = 0 is 


i"{cp 0 ) .o 

o = %o) + " <?cP + ~ kp " ^(P + 


(B3) 


Substituting equation (B2) into (B3) gives 


-% 0 ) f" (<Pq) 

0 = % Q ) + f («p 0 ) (1 + 5) + 

° ° f’(cPn) 2 


'%n)' 


n2 


^ W n) 


(1 + + 


(B4) 


n 

For 1 6 1 « 1, 6* and higher order terms may be neglected, and thus 


6 = 


"<?o>[M ! 


2 {% 0 >[ f > 0>] 2 • f "^ 0 >[ f M 2 } 


(B5) 


Substituting into equation (B2) results in 
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CM 

'"o 


] 2 - i"(<P Q )f(<PQ) 

2 j 

[ f, W 

f - f”(<p 0 )% 0 )j 


(B6) 


which is the Newton-Raphson expression modified. 

For small values of f ' , the Newton-Raphson technique provides slow convergence. 

In such cases, the regula-falsi method was used. The latter technique works best for 
small values of f r . 

In the computer program, 25 iterations were allowed to provide convergence such 

-3 

that the cp calculated in successive iterations differed by less than 10 . A sampling of 

test cases showed that, on the average, three iterations and seldom over five iterations 
were required for convergence when the modified Newton-Raphson technique was used. 

It was rarely necessary to use the auxiliary regula-falsi method for determining <p. 
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